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ABSTRACT 

We report the detection of the double neutron star binary (DNSB) B1534+12 (= J1537-I-1155) with 
the Chandra X-ray Observatory. This DNSB (Porb — 10.1 hr) consists of the miUisecond (recycled) 
pulsar J1537-I-1155A {Pa = 37.9 ms) and a neutron star not detected in the radio. After the remark- 
able double pulsar binary J0737— 3039, it is the only other DNSB detected in X-rays. We measured 
the flux of (2.2 ± 0.6) x 10^^'^ ergs s~^ cm~^ in the 0.3-6 keV band. The small number of collected 
counts allows only crude estimates of spectral parameters. The power-law fit yields the photon index 
F = 3.2 ± 0.5 and the unabsorbed 0.2-10 keV luminosity Lx « 6 x lO^^ergss-^ « 3 x 10"''^;^, 
where Ea is the spin-down power of J1537-f-1155A. Alternatively, the spectrum can be fitted by a 
blackbody model with T « 2.2 MK and the projected emitting area of ~ 5 x 10^ m^. The distribution 
of photon arrival times over binary orbital phase shows a deficit of X-ray emission around apastron, 
which suggests that the emission is caused by interaction of the relativistic wind from J1537-I-1155A 
with its neutron star companion. We also reanalyzed the Chandra and XMM-Newton observations of 
J0737— 3039 and found that its X-ray spectrum is similar to the spectrum of B1534-I-12, and its X-ray 
luminosity is about the same fraction of Ea, which suggests similar X-ray emission mechanisms. How- 
ever, the X-ray emission from J0737— 3039 does not show orbital phase dependence. This difference 
can be explained by the smaller eccentricity of J0737— 3039 or a smaller misalignment between the 
equatorial plane of the millisecond pulsar and the orbital plane of the binary. 

Subject headings: pulsars: individual (PSR B1534-hl2 = PSR J1537-H1155, PSR J0737-3039A, PSR 
J0737— 3039B) — stars: neutron — X-rays: stars 



1. INTRODUCTION 

The sample of known double neutron star binaries 
(DNSBs) has nearly doubled over the last six years. Cur- 
rently, eight such binaries are known (they are shown in 
the P-P diagram in Fig. 1). All of them are found via 
radio timing of pulsars which provides binary orbital pa- 
rameters, including the companion masses. DNSBs are 
thought to be formed in consequent supernova explosions 
in the course of evolution of massive binary systems. The 
details of binary evolution depend on initial star masses 
and binary orbital parameters (see Stairs 2004 for a re- 
view) . It is generally believed that once the first neutron 
star (NS) is formed in the binary, accretion of matter 
from the evolved second star spins up the NS so that 
it becomes a millisecond (recycled) pulsar (MSP). Even- 
tually, the second star explodes as a supernova and, if 
the explosion does not disrupt the binary or disintegrate 
the second star, a DNSB (or a NS-black hole binary) is 
formed. 

Radio observations of DNSBs allow one to measure 
the masses of both NSs and test the predictions of Gen- 
eral Relativity (GR). The recent discovery of the tight 
DNSB J0737-3039 (Porb = 2.45 hours; see Burgay et al. 
2003; Lyne et al. 2004; and Tables 1 and 2 for other de- 
tails), where both NSs are radio pulsars orbiting in the 
plane seen nearly edge-on (the inclination angle i differs 
from 90° by 0?29 ± 0?14; Coles et al. 2005), dramati- 
cally boosted interest in DNSBs. Apart from measuring 
the GR parameters, the radio observations have shown 
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variations of the pulsed flux with orbital phase, for both 
pulsars. The radio pulse of the MSP J0737-3039A (here- 
after J0737A) disappears around its superior conjunc- 
tion^, which has been interpreted as an eclipse of J0737A 
( Lyne et al. 2004), caused by synchrotron absorption in 
the magnetosphere of J0737B (Rafikov & Goldreich 2005; 
Lyutikov & Thompson 2005). Even more interesting is 
the finding that the generally faint pulsed radio emission 
of J0737B is strongly enhanced during two short orbital 
phase intervals of ^ 0.1 of the period around the phase 
of inferior conjunction of J0737B. These changes suggest 
that the magnetosphere of J0737B responds to the wind 
or radiation of the more energetic J0737A (Lyne et al. 
2004; Jenet & Ransom 2004; Demorest at al. 2004; Zhang 
& Loeb 2004). 

The J0737 was the first and, until this work, the 
only DNSB detected in X-rays. It has been observed 
with both Chandra (McLaughhn et al. 2004) and XMM- 
Newton (Campana, Possenti & Burgay 2004; Pellizzoni 
et al. 2004). Its luminosity is Lx ~ (2-3) x lO^" ergs s'^ 
in the 0.2-10 keV band, assuming a 500 pc distance. The 
X-ray spectrum is rather soft; it can be fitted by either 
a power-law (PL) model with a photon index F ~ 3- 
4 or a blackbody (BB) model with kT w 0.2 keV and 
projected emitting area A ~ 3 x 10"* m^. No X-ray vari- 
ability has been reported for J0737 so far. A number 
of different interpretations of the X-ray emission have 
been suggested. These include emission from the sur- 
face and/or magnetosphere of J0737A (similar to solitary 

^ At the superior (inferior) conjunction of J0737A, it is located 
behind (in front of) J0737B along the observer's line-of-sight. 



MSPs), emission from the J0737B's magnetosphere (sur- 
face) energized (heated) by the J0737A's wind (Zhang & 
Loeb 2004), and emission by particles accelerated in a 
bow shock produced by the J0737A's wind in the inter- 
stellar medium (ISM) or in the magnetosphere of JG737B 
(Granot & Meszaros 2004). However, because of the 
small number of X-ray counts collected, it is difficult to 
discriminate between these models from observations of 
just one DNSB. 

In this paper we report the detection of X-rays from 
the second DNSB, B1534+12 (=J1537+1155; hereafter 
J1537), discovered by Wolszczan (1991). One of the two 
NSs in this DNSB is a millisecond pulsar (J1537A here- 
after; period Pa = 37.9 ms, characteristic age ta = 
PaI'^Pa = 248 Myr, spin-down power Ea = 1.8 x 10^^ 
ergs s^^) on a 10.1 hour eccentric (e = 0.274) orbit. As 
one can see from Figure 1 and Table 2, the J1537A's pa- 
rameters are generally close to those of J0737A, although 
J1537A is a factor of 3 less energetic. The dispersion 
measure of J1537A suggests a distance of about 0.7 kpc. 
The comparison of the observed orbital period decay with 
that predicted by GR allowed Stairs et al. (2002) to eval- 
uate the distance more accurately: P) = 1.02 ± 0.05 kpc. 
Multi-year timing observations in the radio (Stairs et al. 
2004) suggest that we are looking at J1537 close to its 
orbital plane (i ~ and J1537A is an almost or- 

thogonal rotator (i.e., its magnetic axis and the line of 
sight are nearly orthogonal to the spin axis: a « 103°, 
^ « 97° at the epoch of our observation^) whose spin 
axis is misaligned with the orbital angular momentum 
by (5 = 25° ±4° (Stairs et al. 2004; Thorsett et al. 2005). 
The second NS (hereafter J1537B) has not been detected 
as a radio pulsar, but its mass, Mb — (1.345±0.001)Mq, 
has been measured (Stairs et al. 2002). The detection of 
X-rays from J0737 prompted us to observe J1537, the 
second nearest DNSB. We describe the Chandra obser- 
vation of J1537 in §2 and the observational results in §3, 
which also includes our reanalysis of the X-ray observa- 
tions of J0737 and a comparison of the X-ray properties 
of these DNSBs. We discuss plausible X-ray emission 
models and compare their predictions with the observed 
X-ray properties in §4 and summarize our results in §5. 

2. OBSERVATION 

We observed J1537 with the Advanced CCD Imaging 
Spectrometer (ACIS) aboard Chandra on 2005 April 10 
(start time 53,470.352222 MJD). The useful scientific ex- 
posure time was 36,080 s (first 935 s of the 37,015 s to- 
tal exposure are excluded from the Good Time Intervals 
[GTIs]). The observation was carried out in Very Faint 
mode, and the pulsar was imaged on S3 chip of the ACIS- 
S array. The detector was operated in Full Frame mode 
which provides time resolution of 3.2 seconds. The data 
were reduced using the Chandra Interactive Analysis of 
Observations (CIAO) software (ver. 3.2.1; CALDB ver. 
3.0.3). 

3. X-RAY IMAGE AND SPECTRUM 

Figure 2 shows the ACIS-S3 image of the J1537's field. 
The X-ray source is clearly seen at R.A. = 15'^37'"09!986, 
Decl. = -f 11°55'55'.'39 (the combined aspect determina- 
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Fig. 1. — P-P diagram for Ri 1400 radio pulsars (small dots) 
from the ATNF catalog (Manchester et al. 2005). Lines of constant 
pulsar age, magnetic field, and E are shown. Eight DNSBs are 
marked by different symbols. J1537A is marked by the filled star. 
The J0737 components are marked by open stars connected with 
the straight dashed line. The rest of the DNSBs are marked by 
filled circles. See the electronic edition of the Journal for the color 
version of this figure. 
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Fig. 2.— ACIS-S3 image of the J1537 field. The circle of 1'.'23 
radius shows the radio position of PSR J1537-I-1155 from Stairs et 
al. (2002). 

tion and centroiding uncertainty is 0'.'7 at a 90% con- 
fidence level). Since this position is within 0'.'5 of the 
J1537A's radio position (Stairs et al. 2002), we conclude 
with confidence that we detected the X-ray emission from 
J1537. The distribution of source counts in the ACIS im- 
age is consistent with that of a point source. 

We extracted the pulsar spectrum from a circular aper- 
ture of 2.5 ACIS pixels radius (= 1'.'23; « 90% encircled 
energy radius) using the CIAO psextract task. The 
background was extracted from a 5" < r < 22" annu- 
lus centered on the source. The total number of counts 
within the source aperture is 16, of which 98.4% are ex- 
pected to come from the source. The observed source 
fiux is (2.2 ± 0.6) X 10-1^ erg s^^ cm'^ in the 0.3-6 
keV range. We group the counts into 4 spectral bins, 
with each bin having 3-5 counts (see Fig. 3), and fit the 
spectrum with absorbed BB and PL models using the C- 
statistic (Cash 1979) implemented in XSPEC, ver. 11.3.0 
(the fitting parameters are given in Table 3). To obtain 
constrained fits, we had to freeze the hydrogen column 
density, TVh- The pulsar's dispersion measure, DM= 11.6 
cm~^ pc, corresponds to iVn = 3.6 x 10^° cm~^ (assum- 
ing a 10% ISM ionization). This value is close to the 
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Fig. 3. — X-ray spectrum of J1537 fitted with PL (top) and 
BB {bottom) models. The contributions of the energy bins into 
the best-fit C-statistic are shown in lower panels, multiplied by —1 
when the number of data counts is smaller than the number of 

model counts. 

total Galactic absorption in neutral hydrogen, 3.6 x 10^'^ 
cm^^ (Dickey & Lockman 1990) in the direction of J1537 
(/ = 19?85, b = 48?34), which is located well above the 
Galactic plane {z = 0.76 kpc at the distance of 1 kpc). 
Therefore, we adopt A^h = 3.6 x 10^° cm~^ in the X-ray 
fits. 

The spectrum fits a soft PL with a photon index 
r = 2.7-3.7 and an isotropic unabsorbed 0.2-10 keV 
luminosity Lx = ^ttD'^F^''^^^ ~ i^-^tH) x lO^^ ergs 
s^^ at D = 1 kpc. The temperature and the pro- 
jected area of the emitting region obtained from the BB 
fit are strongly correlated (see Fig. 4), which results in 
large uncertainty of these parameters. The BB model 
gives the best- fit temperature Tbb = 2.2 MK and a pro- 
jected emitting area ^ = 5 x 10'^ m^. The projected 
area is much smaller than that of a NS, ~ 3 x 10® m^, 
suggesting that the thermal radiation might be emit- 
ted from small, hot polar caps with an effective ra- 
dius R ^ 39(cos<;)~^/^ m and a bolometric luminos- 
ity Lboi = aT^A{cos^)-'^ = (6.8 ± 2.1) x 1028(cos?)-i 
ergs s~^, where (cos?) is the time-averaged cosine of the 
angle between the line of sight and the normal to the 
visible polar cap ((coscr) ~ 0.62 for the axis orienta- 
tions suggested by Stairs et al. 2004, assuming a cen- 
tered dipole with two equivalent polar caps and neglect- 
ing the GR bending of photon trajectories). The X-ray 
luminosity of J1537 corresponds to 3.4 x 10~^£'^ and 
3.8 X 10-5(cos<?)-1£a = 6.1 X lO'^E for the PL (0.2-10 



keV band) and BB fits, respectively. The scarce statis- 
tics does not allow us to discriminate between the BB 
and PL models. 

To compare the spectra and light curves of J 1537 and 
J0737, we reanalyzed the Chandra ACIS and XMM- 
Newton EPIC observations of J0737, using the same 
spectral models and the same approach to the data anal- 
ysis as for J1537. The observations have been described 
by McLaughhn et al. (2004), Campana et al. (2004), and 
Pellizzoni et al. (2004). The results of our fits are pre- 
sented in Table 3 and Figure 5. For the PL fit of the 
ACIS spectrum, we obtained about the same parame- 
ters as McLaughlin et al. (these authors did not fit the 
BB model). As for the EPIC data, we found, in ac- 
cordance with Campana et al., that the EPIC-pn data 
are too noisy for a useful analysis because of the very 
high background in the timing mode. Therefore, we only 
used the MOSl and M0S2 data. For each of the MO'S 
detectors, we created a light curve with binning time 
of 10 s and constructed GTIs excluding the time inter- 
vals during which the total chip count rate exceeded 3 
counts s^^ (the remaining useful time is 24.62 and 24.60 
ks for MOSl and M0S2, respectively). The spectra we 
extracted from circular apertures of 46" radius (^ 90% 
encircled energy). Comparing the properties of the X-ray 
radiation from J0737 and J1537, we see that the spectral 
parameters are very similar: F ss 3 for the PL fit, and 
kT 0.2 keV for the BB fit. 

In the upper panel of Figure 6 we show the distribu- 
tion of arrival times for the 16 detected photons over 
binary orbital phases of J1537 {(poih = corresponds 
to periastron passage)^. Notice that no events arrived 
during the 12,795 s interval between MJD 53470.46046 
and 53470.60855, which corresponds to the phase inter- 
val [0.365; 0.717] (A(?!)orb = 0.352). This interval includes 
the phase of apastron, (f) = 0.5, and the phase of su- 
perior conjunction of J0537A, (f) = 0.408 at the epoch 
of our Chandra observation. The probability of getting 
zero counts by chance in the 12,795 s interval is 0.0034 
(assuming Poisson statistics), i.e. the observed "gap" has 
a w 3 cr significance. For comparison, we plotted similar 
distributions for the Chandra and XMM-Newton obser- 
vations of J0737 (middle and bottom panels of Fig. 6), 
which do not show a statistically significant dependence 
on orbital phase (Fig. 7). 

4. DISCUSSION. 

The small number of detected counts does not allow 
one to determine the nature of the X-ray emission from 
J1537 unequivocally. Among viable models are intrinsic 
emission from the surface or magnetosphere of the NS(s) 
(just as in solitary pulsars), the synchrotron emission 
resulting from the interaction between the pulsar winds 
or between the wind of J1537A and the magnetosphere of 
J1537B, and emission from the inner magnetosphere or 
surface of J1537B triggered by the relativistic particles 
from the J1537A's wind. Below we briefly discuss these 
interpretations, with emphasis on different dependences 
of the X-ray flux on binary phase. 

To calculate the event orbital phases, we used the binary 
cphcmcris for the epoch of our observation kindly provided by 
Ingrid Stairs; Porb = 0.420737299122 days, epoch of periastron 
53470.309046790 MJD, longitude of periastron 290?004585816. 
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Fig. 4.— Confidence contours (68% and 90%) for the PL (top) 
and BB (bottom) fits to the spectrum of J1537. The PL normal- 
ization is in units of 10~^ photons cm~^ s~^ keV~^ at 1 keV. The 
BB normalization (vertical axis) is the projected emitting area in 
units of m^, for D = 1 kpc. The lines of constant unabsorbed flux 
(top panel; in units of 10~^^ ergs cm"'^ s~^, in 0.2—10 keV band) 
and constant Li,oi(cos'j) (bottom panel; in units of 10^* ergs s~^) 
are plotted as dashed lines, for fixed Nn = 3.6 X 10^'' cm~^. 



4.1. Phase- dependent X-ray emission powered by the 
wind from pulsar A 

The apparent lack of J1537's X-ray emission near the 
apastron suggests that a phase-dependent interaction 
mechanism is responsible for the X-ray emission in close 
DNSBs. The interaction can proceed via several different 
mechanisms, and the pulsar winds are expected to play 
an important role in each of them. Below we consider 
plausible scenarios. 

4.1.1. Interaction of pulsar A's wind with the outer 
magnetosphere of pulsar B 

As discussed in many papers on pulsars in bina- 
ries (e.g., Arons & Tavani 1993), phase-dependent syn- 
chrotron emission can be produced when the pulsar's 
wind collides with the wind or magnetosphere of the sec- 
ondary companion. For a DNSB, the position of the 
interaction site (e.g., of the bow-shock head) on the line 
connecting the two pulsars can be estimated by balanc- 
ing the pressure of the pulsar A's wind with the pressure 



Fig. 5. — Same as in Fig. 4 but for J0737. Solid and dash- 
dotted contours are from the XMM-Newton EPIC MOS and Chan- 
dra ACIS observations, respectively. See the electronic edition of 
the Journal for the color version of this figure. 

of pulsar B's wind or magnetosphere: 



EAfA 



Ef 



47rc(d 
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rB)^c A-Kcrj. 



' fs if TB > Rlc.B 

93 (^)'if TB-^RiJ^^ 



where dAB is the distance between the binary compan- 
ions, rB is the distance between the interaction site and 
pulsar B, and Ric,b — 4.8 x IO^Pb cm is the light cyhn- 
der radius for pulsar B (approximately equal to the size 
of undisturbed magnetosphere). The factors and 
/s take into account possible anisotropy of the pulsar 
winds (/ = 1 for an isotropic wind). They depend on 
the colatitudes, and of the wind interaction site 
with respect to the A's and B's spin axes because the 
pulsar wind should be axially symmetric; for instance, 
/ = (3/2)sin^t9 for an angular distribution similar to 
that of magnetic dipole radiation (cf. Arons et al. 2004), 



/ 



(sinA)-i if |i? - 7r/2| < A 
otherwise 



(2) 



for an equatorial outflow with a uniform energy flux 
within sharp boundaries. The factor gB — (3/4) (1 + 
3 cos^ Qb) takes into account the dependence of the mag- 
netic pressure on the polar angle 8 b with respect to the 
B's magnetic axis. 
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Fic;. 6. — Distributions of photon arrival times over binary 
phase for the Chandra ACIS observations of J1537 (top) and J0737 
(middle), and for the XMM-Ncwton MOSl and MOS2 observation 
of J0737 (bottom). Zero phase corresponds to periastron passage. 
Only first 14.05 ks of the XMM-Ncwton exposure, least affected 
by the flare background, are used here. The hatched area in the 
top panel shows the time interval excluded from the GTIs. The 
dashed vertical lines show the phases of superior (S) and inferior (I) 
conjunctions for J1537A and J0737A at the epochs of the Chandra 
and XMM-Newton observations. 




Fig. 7. — X-ray light curve of J0737 folded with the binary 
period. We used both Chandra ACIS and XMM-Ncwton MOS 
data taken over the length of a single binary revolution in each case 
(most of the remaining 51-ks-long XMM-Newton exposure suffered 
from strong background flares.). The horizontal dotted line shows 
the background level. 
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Fig. 8. — Configuration of the J1537 binary at the epoch of 
the Chandra observation (53,470 MJD). View from above the or- 
bital plane with the Earth to the right is shown. The ellipses show 
the orbits of J1537A and J1537B around the center of mass of 
the system. The filled and open circles show the positions of the 
components at superior and inferior conjunction of J1537A. The 
vertical solid line is the line of nodes, and the dashed straight line 
is the apsidal line. The ticks on the J1537B's orbit show the distri- 
bution of the 16 detected photons over the orbit, while the shaded 
part of the orbit (around the line of nodes) corresponds to the ex- 
cluded time interval at the beginning of the observation. (Notice 
that the distribution of photons over the orbit looks different from 
that over phases in Fig. 6 because of the substantial eccentricity; 
in particular, the gap near the periastron is much larger because of 
the increased orbital velocity.) The shaded sectors show the ranges 
of true anomalies (counted clockwise from periastron) in which 
the lack of detected photons could be explained by a misalignment 
of the A's spin axis and the binary angular momentum, assuming 
an equatorial outflow of the A's wind (see §4.1.1). 

Since B is expected to be a slowly rotating, old (per- 
haps even dead) pulsar, we can assume Eb <^ Ea (e.g., 
Eb/Ea = 2.8 X 10"^ for J0737), which gives 

f dAB{EBfB/EAfAf'^ for TB > Rlc,B , . 

~ I d]l^R^/%iEBgB/EAfA)'/^ for rs < Ric,b ^ ' 

The interaction may occur outside or inside the B's light 

cylinder, dc^pending on the separation and the unknown 
properties of pulsar B. In particular, rs < Ric,b (i-e., 
the B's magnetosphere is compressed by the A's wind on 
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the "daysidc" ) if 

Pb > 0.66 dABMEB,30gB/EA,33fAy^^ s, (4) 
which corresponds to 

r0.3(W/A)'/'s forJ0737 , . 

^ \ 1.1 {dAB/a){EB,sogB/fA)'/^ s for J1537 ' 

where a is the scmimajor axis of the binary relative orbit 
(a = 2.28 X 10" cm and 8.78 x lO^o cm for J1537 and 
J0737, respectively). 

The shocked A's wind can generate synchrotron X-ray 
emission if the wind particles are sufficiently energetic 
and the magnetic field B is strong enough in the emission 
region: 

7^5 > 10^2 (£/5keV), (6) 

where 7 is the Lorentz factor, and E is the photon en- 
ergy. For the plausible case rs < Ric,b, the magnetic 
field at the interaction site can be estimated from the 
pressure balance between the A's wind and the B's mag- 
netic pressure: 

/ • \ 1/2 

/2^a/a\ __a_A/2(7.1G forJ0737,. 
V y "^AB^^ 11-5G forJ1537-^^^ 

At such fields, we need electrons with Lorentz factors 
> 10^-10^ to produce synchrotron X-rays. The recent 
pulsar wind models (e.g., Kirk & Skjairaascn 2003 and 
references therein) suggest that the initially Poynting- 
flux-dominated wind of pulsar A converts most of its 
electromagnetic energy into particle kinetic energy at a 
distance much greater than the binary separation, which 
implies that the bulk Lorentz factor upstream of the in- 
teraction site remains the same as acquired in the pulsar 
A's magnetosphere. Although the predictions of the cur- 
rent pulsar models are rather uncertain in this regard (in 
particular, the Lorentz factor depends on poorly known 
"multiplicity" of the pair cascades; e.g., Hibschman & 
Arons 2001), the required values of 7 are likely too large 
to be produced in the A's magnetosphere, which means 
that additional acceleration must occur at the interac- 
tion site (e.g., caused by reconnection of the magnetic 
fields of the A's wind and B's magnetosphere, similar to 
the magnetic reconnection at the interaction of the solar 
wind with the Earth's magnetosphere). 

If the shocked/accelerated A's wind is the sole source 
of the X-ray emission, the X-ray luminosity should be a 
fraction of the A's wind power: 

Lx = eneeEradex^^A (8) 

where eo ~ fA{r±/'2dAB)'^ is the fraction of A's power 
intercepted by the B's magnetosphere with a transverse 
size r±_ that depends on geometry of the interaction re- 
gion (crudely, rx ~ Ric,b), £e is the fraction of energy of 
the shocked/ accelerated wind carried by particles (elec- 
trons and/or positrons), erad ~ niin(l; fflow/irad) is the 
radiative efficiency, i.e. the fraction of intercepted en- 
ergy emitted as (synchrotron) radiation (frad and tfiow 
are the time of synchrotron losses and the flow time in 
the emission region), and ex is the fraction of synchrotron 
power emitted in the X-ray range. The intercepted frac- 
tion can be estimated as en ~ 6 x 10~^ f Air ±/Ric,b)^ 
and ~ 5 X 10-4/A(r_L/10^°cm)2(a/rfAB)^ for J0737 and 



J1537, respectively. The fraction eg depends on accel- 
eration mechanism (e.g., magnetic reconnection); the 
latter is currently unknown, but the assumption that 
ee is not much lower than unity seems plausible. The 
fraction erad depends on the particle energy spectrum, 
geometry of emitting region, and flow velocity. For a 
particle of given energy, the time of synchrotron losses 
is Trad = 5107^^i?~^s, while the flow time is tnov, ~ 
V^flow ~ 0.3Zio/3~^ s, where / = lO^^/iocm is the travel 
length, and /3 = Vfio^/c. This gives erad ^ min(l;7 x 
10~^liojeB^(3^^). Wc see that electrons emit a very 
small fraction of their energy unless the travel length 
substantially exceeds a typical size of the B's magneto- 
sphere. Finally, the fraction ex cannot be reliably esti- 
mated without knowing the particle energy distribution. 
However, even for ex (and eg) close to unity, the above 
estimates of eq and erad suggest the A's wind should 
be substantially anisotropic (/a ^ 1) and/or the emit- 
ting particles should spend a long time (> 10^ s) in the 
emitting region to explain the ratio Lx/Ea ^ 3 x 10~^ 
observed^ in J0737 and J1537. 

The interaction mechanism and the properties of the 
A's wind can be inferred from the analysis of the de- 
pendence of the X-ray emission on orbital phase. The 
collision of the A's wind with the B's magnetosphere is 
expected to form a bow-shock, convex toward pulsar A, 
with the shocked A's wind flowing in a magnetosheath 
engulfing the B's magnetosphere (Arons et al. 2004; Lyu- 
tikov 2004). It is tempting to explain the X-rays observed 
from J0737 and J 1 537 as synchrotron emission from the 
shocked A's wind in the magnetosheath (as discussed by 
Granot & Meszaros 2004). Such an explanation implies 
an orbital dependence of the X-ray flux caused by the 
relativistic beaming (Doppler boost), showing a mini- 
mum flux at the phase of inferior conjunction of pulsar 
A, when the wind in the sheath is flowing away from 
the observer, and a maximum flux around the superior 
conjunction. For a reasonable flow speed, vnow = l3c, 
the Doppler boost modulation can be quite large, up to 
i^max/-Fmi„ = [(1 + /3)/(l - /9)]^+^ especially at the ob- 
served large photon index, F ~ 3. However, the observed 
distributions of photons over orbital phase do not show a 
statistically significant flux deflcit at the phase of A's in- 
ferior conjunction, in neither J0737 nor J1537 (see Figs. 6 
and 7). Therefore, such an explanation can be ruled out. 
It seems surprising that no orbital modulation caused by 
the anisotropy of the flow in the putative magnetosheath 
is seen (unless, of course, our assumption, that the X-rays 
are produced by interaction of the A's wind with the B's 
magnetosphere, is wrong). We can only speculate that 
either the magnetosheath does not form because the A's 
wind is strongly magnetized or the X-ray emitting par- 
ticles accelerated by magnetic reconnection do not flow 
into the magnetosheath. 

Having ruled out the magnetosheath as the source of 
the observed X-rays, we can, nevertheless, suggest two 
explanations for the observed orbital dependence of X- 
ray emission in J1537, such that they do not contradict 
to the lack of such dependence in J0737. First, this dif- 
ference may be connected with the fact that J1537 has a 

We should note that the Lx estimates in §3 are based on the 
assumption of isotropic X-ray emission, which is not necessarily 
true; the actual Lx can be lower or higher than estimated. 
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substantial eccentricity while J0737 is on an almost circu- 
lar orbit. At apastron, the separation cIab is larger than 
at periastron by a factor (H-e)/(l — e), which is 1.755 for 
J1537, vs. 1.192 for J0737. The decrease of the magnetic 
field and particle number density, rig, at the interaction 
site with increasing separation should lead to a lower X- 
ray luminosity at apastron in eccentric DNSBs. Crudely, 
assuming ex does not depend on (Iab substantially, the 
luminosity behaves as Lx ~ ^xneVVsyn oc Vd^g, where 
Vsyn ex oc d~^g is the total synchrotron power per elec- 
tron , and V is the emission volume. The dependence of 
the latter on dAB is rather uncertain, but it seems reason- 
able to assume V (X r^. Under this assumption, we ob- 
tain V oc d^g or y oc dAB for vb > Ric,B or tb < Ric,B, 
respectively, so the increase of V with dAB does not com- 
pensate for the decrease caused by decreasing B and rig. 
In particular, Lx oc d'^g for tb < Ric,B, which means 
that we can expect the luminosity at apastron to be a fac- 
tor of 5.4 lower than at periastron, for J1537. Although 
the effect of varying binary separation can explain the 
observed deficit around apastron in J1537, we do not see 
an increased count rate at periastron, predicted by this 
model; however, we cannot rule it out because of the 
scarce statistics of our data. 

In addition to the varying binary separation, the phase 
dependence of the X-ray flux in J1537 can be caused by 
a misalignment between the binary orbital plane and the 
equatorial plane of pulsar A, where most of the A's wind 
is expected to be confined. In this case, one should ex- 
pect brighter X-ray emission at orbital phases when pul- 
sar B plunges in the dense wind flowing in the equatorial 
plane of pulsar A. For instance, if the equatorial outflow 
is conflned between two conical surfaces, within an angle 
±A from the equatorial plane, and the A's spin axis is 
inclined at an angle 5 to the direction of the orbital angu- 
lar momentum, then, at (5 > A, the wind "misses" pulsar 
B in two segments of the orbit: < r; and \tv — < t/, 
where tp is the azimuthal angle counted from the projec- 
tion of the A's spin onto the orbital plane, and ry is given 
by the equation 

cot?7 = tanA (tan^ 5 - tan^ A) "^/^ (cos 5)"^ . (9) 

In other words, one should expect a photon deficit in two 
phase intervals corresponding to two sectors symmetric 
with respect to the A's spin projection onto the orbital 
plane. If such a phase dependence is observed, then rj 
can be measured, and A can be found from the equation 

tan^ A = cot^ 77 sin^ 5{l + cot^ ?? cos^ 5)"^ (10) 

For J1537, we determined the true anomaly for each 
of the 16 events and marked the events on the J1537B 

orbit in Figure 8. Because of the higher (lower) NS ve- 
locities near the periastron (apastron), a relatively small 
phase segment around (/> = with no photons detected, 
—0.114 < (j) < 0.043, translates into a sector of true 
anomalies, -68° <e < 28°, whose width, A6'per = 95°, 
is even larger than that of the photon-free sector at apas- 
tron, A6lap = 80° (150° < e < 231°). If we inter- 
pret the photon deficit in these sectors as causcxl by the 
anisotropy of the wind outflow, then, from the symme- 
try requirement, the sectors in which the J1537A's wind 
misses J1537B (shaded in Fig. 8) are -30° < 6I < 28° 
and 150° < 6 < 208°, i.e., the projection of the A's spin 
onto the orbital plane almost coincides with the apsidal 



fine, and r/ < 29°. Since S = 25° ± 4° for .11537 (Stairs 
et al. 2002), equation (10) gives a lower limit on the 
half-opening angle of the equatorial outflow, A > 18°- 
25°. Being consistent with the observed phase depen- 
dence, such an explanation is also attractive because it 
implies an increased particle and energy supply in the 
two segments of orbit where the J1537A's wind encoun- 
ters J1537B [Ja can be as large as 2.4-3.2, according to 
eq. [2]), a higher magnetic field at the interaction region 

(by a factor of /^^^), and, correspondingly, a smaller ra- 
tio tfiow/^rad (by a factor of /a)- These consequences of 
the wind anisotropy alleviate the energetic and geomet- 
ric restrictions inherent in the interpretation of the X-ray 
emission as caused by the interaction of the A's wind with 
the B's magnetosphere. We should emphasize, however, 
that the X-ray flux deficit around periastron is not statis- 
tically significant, and this hypothesis can be confirmed 
only by a deeper observation. For J0737, the angle S is 
rather uncertain (0 < ^ < 60°; Manchester et al. 2005), 
so one can assume that the lack of a clear orbital phase 
dependence means that 6 < A, i.e., the J0737A's wind 
interacts with J0737B throughout the whole orbit. 

In addition to the energetics and geometry (luminos- 
ity and its phase dependence), some information on the 
interaction mechanism can be obtained from the X-ray 
spectra. In this regard, we should note that the spectra 
observed in J1537 and J0737 are surprisingly soft (F ^ 3) 
in comparison with, e.g., the spectra observed in pulsar 
wind nebulae (F ~ 1-2). The models of particle ac- 
celeration by ultrarelativistic MHD shocks also suggest 
much flatter spectra, F ~ 1.6, corresponding to the slope 
p ~ 2.2 of the particle energy spectrum (Achterberg et al. 
2001). We can speculate, however, that acceleration via 
relativistic magnetic reconnection might produce softer 
spectra. Also, we do not exclude the possibility that the 
observed spectra consist of two components: a hard PL 
component from the shocked wind and a soft thermal 
component from the polar caps of pulsar A (see §4.2). 

4.1.2. Emission from B's inner magnetosphere and/or 
surface induced by captured A 's wind 

There is a possibility that some fraction of the A's 
wind is captured by the B's magnetosphere (e.g., simi- 
lar to capturing solar wind in the Earth magnetosphere). 
The captured wind can increase the magnetospheric pair 
density or even precipitate onto the NS surface (such a 
possibility has been discussed by Zhang & Loeb 2004 to 
explain the flaring radio emission of J0737B). The cap- 
tured relativistic particles lose a substantial fraction of 
their energy to synchrotron and inverse Compton radia- 
tion, thereby increasing the radiative efficiency, erad, of 
the intercepted wind. The particles traveling along the 
open fleld lines may reach the NS polar caps and heat 
them to X-ray temperatures, up to a few million kelvins, 
so that we could see not only nonthermal but also ther- 
mal X-ray emission, which might explain the softness of 
the X-ray spectra observed in J0737 and J1537. The 
total luminosity from these emission mechanisms is, of 
course, a small fraction of Ea', it can be described by an 
equation similar to equation (8) but with an additional 
factor in its right-hand side, Ccapj which is the fraction 
of intercepted wind particles captured by the B's magne- 
tosphere. Similar to the emission produced by the wind 
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interaction with the outer magnetosphcrc (§4.1.1), the X- 
ray luminosity should decrease with increasing cIab (i-e., 
it should be lower near apastron) and vanish during the 
phase intervals when pulsar B is far from the equatorial 
plane of pulsar A^. In addition, the fraction Ccap should 
depend on the angle between the B's magnetic axis and 
the line connecting the two pulsars. This may lead to 
an additional modulation of Lx with the orbital period 
and the period of pulsar B unless the B's magnetic axis 
is perpendicular to the orbital plane, as expected if the 
main braking mechanism for pulsar B is its interaction 
with the A's wind (Arons ct al. 2004; Lyutikov 2004). 
In addition, the X-ray emission caused by the captured 
wind should show pulsations with the B's spin period. 

There is also an interesting possibility that the wind's 
particles that reach the B's surface and heat the polar 
caps can facilitate the extraction of electrons or ions from 
the polar cap surface by the pulsar's electric field, trig- 
gering the pair cascade in pulsar B. In this case, the A's 
wind works as a catalyst of B's X-ray (and possibly radio) 
emission, providing an additional luminosity associated 
with the loss of the pulsar B's spin energy^. The orbital 
phase dependence of this luminosity should be similar to 
that of the above-discussed captured wind. In particu- 
lar, even if B is a dead pulsar when no wind blows on 
it (e.g., when it is far from the A's equatorial plane), it 
may become "resurrected" on the orbit segments where 
it captures the wind. 

Quantitative predictions of the properties of the X-ray 
emission induced by the captured wind are highly uncer- 
tain, depending on the unknown wind properties, struc- 
ture of the magnetic field near the interaction region, etc. 
We can crudely estimate the number of A's wind parti- 
cles captured by B per unit time: SNa = K^A^GJ.ACnecap; 
where TVqj.a ^ {3EAc/2e'^)'^/'^ is the Goldrcich- 
Julian rate of particle ejection from the A's magneto- 
sphere, KA the pair production multiplicity, and eo ~ 
fAirs/'i'dAB)^- Assuming tb < Ric.B, we obtain SNa ^ 
2.2 X WKAecaAa/dAB)^/'PT{fA9BEB^3oy^'s-^ for 
J1537, and 5Na ~ 4.4 x 102^K^ecap(/l5B)^/^ s'^ for 
J0737. If the captured wind is the only source of X- 
rays, the product KA^cap should be large enough to 
provide the observed X-ray luminosity: e.g., /CACcap > 

3 X lOHflgB)-'/'Ps'^'EB%\o,/dAB)-'/HlO''m,^cy£) 
in J1537 and KACcap > 0.8 x lO^iflgBy^/'HW^'rneC^/S) 
in J0737, where S is an average energy per captured par- 
ticle rcradiated in X-rays (Zhang & Loeb 2004 assume 
KA ~ 10^, Ccap 0.1 in their estimates for J0737; as we 
have mentioned above, at such large ka no X-rays can 
be expected from the wind that was not captured). For 
Ka ~ 10^ obtained in standard cascade theory, we have 
to invoke the above-mentioned activation of pulsar B by 
the captured wind to explain the observed energetics. 

Note that a large value of KA^cap is also required for 
SNa to exceed the intrinsic B's pair production rate, 

* We should mention that such orbital phase dependencies axe 
expected not only for the nonthermal but also for the thermal com- 
ponent because the polar cap cooling time scales are very short 
(about a few microseconds for B ~ 10^^ G; e.g., Gil et al. 2003). 

® If such a mechanism works, we should expect that Eg (and 
Pb ) are increased during the "activation phases" . 



Nb — kbNqj^b — 4: A X 10'^^ Eg^f^ Kb s ^. It occurs when 

«Aecap/«B > 200 {dAB/ar/-'PB^^'iflgB)-'/'EB%^^' 

and 13 (/Iffs)"^/^ for J1537 and J0737, respectively. 

Thus, we can conclude that energizing the B's mag- 
netosphere or heating its surface by the particles from 

the A's wind can explain the observed dependence of the 
X-ray flux on binary phase. However, it would require a 
rather high pair production multiplicity in pulsar A and 
considerable anisotropy of the A's wind. 

4.2. X-ray emission from individual pulsars 

Since solitary pulsars are known to emit X-rays, at 
least some fraction of the detected X-rays could come 

from one or both of the DNSB components, even if 
the X-ray emission caused by the interaction of the A's 
wind with the B's magnetosphere were negligible. In 
particular, MSPs can exhibit both nonthermal (mag- 
netospheric) and thermal (polar cap) components in 
their X-ray spectra, the latter dominating in MSPs with 
E < 10"^^ ergs s^^. Typical thermal conversion efficien- 
cies, r/th = Lboi/E, are (0.3-3) x 10"^ (Zavlin 2006, 
and references therein). The efficiencies inferred from 
the BB fits of the J1537 and J0737 spectra, 6.1 x lO''^ 
and 4.6 x 10~^(cos<j)~^D|go, respectively, are within 
that range^. The polar cap radii derived from the BB 
fits, are smaller than the conventional PC radii [rpc = 
(27ri?^^/cP)^/^ = 740 m and 960 m, respectively), sim- 
ilar to other pulsars (e.g., Zavin & Pavlov 2004; Za- 
vlin 2006). Thus, the luminosities and the soft spectra 
of J1537 and J0737 are consistent with being emitted 
from polar caps of J1537A and J0737A, but of course 
such emission would not show any dependence on binary 
phase. 

The X-ray emission in J1537 might also originate from 
J1537B if this NS is not a "dead pulsar", i.e. it is still 
capable of producing e+e~ pairs emitting curvature and 
synchrotron radiation in the magnetosphere and heating 
its polar caps. Indeed, the X-ray spectrum of J 1537 re- 
sembles the spectra of old (r ~ 1-10 Myr) solitary radio 
pulsars (e.g., PSR B0950-I-08; Zavhn and Pavlov 2004). 
However, the J1537's luminosity (e.g., Lx = 8.5 x 10^^ 
ergs s^^ in the 1 TO keV band, for the PL fit) is rather 
high for such an old pulsar* ; it would either imply a high 
spin-down luminosity, SJ1537B > 10^^ ergs s~^, or sug- 
gest that the X-ray efficiencies of old pulsars strongly 
grow with decreasing E, which is hardly supported by 
observations (cf. Fig. 5 in Kargaltsev et al. 2006). More- 
over, such an interpretation is firmly excluded for J0737, 
whose X-ray luminosity exceeds i^jo737B = 1.6 x 10^° ergs 
s-i. 

Thus, if the observed J1537's emission originates from 
the surface or magnetosphere of one of the binary com- 
panions, J1537A is a much more viable candidate than 
J1537B. The assumption that most of the observed X- 

ray emission is generated at the polar caps (or magneto- 
sphere) of J1537A would be strongly supported if X-ray 

^ Since the orientations of the axes is poorly known for J0737A 
(Demorcst ct al. 2004; Manchester et al. 2005), the value of (cos?) 
is very uncertain. 

Based on the J1537's proper motion and distance from the 
Galactic plane, Thorsett et al. (2005) argue that the SN explosion 
that produced J1537B occurred > 10 Myrs ago. 
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pulsations at the J1537A period are found in future ob- 
servations. An indirect support for this interpretation 
would also be provided by detection of X-ray pulsations 
from J0737 at the J0737A's period. On the other hand, 
confirming the reduced count rate near apastron would 
virtually rule out this interpretation and put an interest- 
ing upper limit on the MSP luminosity. 

4.3. Emission from the A \s wind interacting with the 

interstellar m,edium 

The velocities of J1537 and J0737 in the plane of 
the sky, v±_ ~ 120 and 70 km s^^, respectively, are 
large enough to assume that these DNSBs are moving 
through the ISM with supersonic speeds, i.e., v > Cs = 
15(/i/0.6)-i/2ry^ km s-\ where /i and T = IQ^T^ K 
are the molecular weight and temperature. At a su- 
personic speed, the pulsar's wind forms a termination 
shock (TS) in the ISM. For an isotropic wind, the TS 
acquires a bullet-like shape (Bucciantini, Amato, & Del 
Zanna 2005), with a distance i?h — {Ea/^t^cpv^Y^'^ = 
4.0 X 10^^£^y33n~^/^wf ^ cm between the pulsar and the 

bullet head [here n = p/(1.66 x 10^^^ g) is the ISM num- 
ber density in atomic mass units per cm^, and v-j is the 
pulsar speed in units of 100 km/s]. For large enough 
in the shocked plasma (see equation [6]), it can generate 
X-rays, as discussed by Granot & Meszaros (2004) for 
J0737. 

The interpretation of the observed X-rays as gener- 
ated by the A's wind interaction with the ISM rather 
than with the B's magnetosphere (§4.1) has two appar- 
ent advantages for explaining the energetics (but not 
the putative dependence on binary phase in J1537): a 
much larger fraction of the A's wind is intercepted (e.g., 
en ~ 0.5), and, possibly, a larger fraction of the A's 
wind energy is carried by particles (cg is closer to unity) 
because a larger fraction of the wind's magnetic energy 
can be converted into particle energy at such large dis- 
tances, i?h ^ dAB- On the other hand, this interpre- 
tation may have its own problems, even with the X-ray 
energetics. First of all, since the wind's magnetic field is 
miich lower at the TS than at the site of interaction with 
the B's magnetosphere, the synchrotron power and radia- 
tive eflaciency, Cradj are much lower, which requires much 
higher particle energies for the radiation to be emitted 
in the X-ray range. Assuming that the magnetization 
parameter (i.e., the ratio of the Poynting flux to the ki- 
netic energy flux) is small, a <C 1, we can, following 
Kennel & Coroniti (1986), estimate the post-shock mag- 
netic field at the shock head: B ^ 3(T^/^(47rpw^)^/^ = 

14 ai^^n^/^i'T /tiG, where ct_2 — cr/lO^^. For the syn- 
chrotron radiation in such a field to be emitted in the 
X-ray range, we should have enough electrons in the en- 
ergy range 7 ~ (1-3) x \Q^aZ.y^'n~^/'^v^^^'^ . Such en- 
ergies are too high to be acquired in the pulsar mag- 
netosphere. The wind electrons can be further accel- 
erated in the pre-shock wind by the conversion of the 
electromagnetic wind's energy into kinetic energy. How- 
ever, irrespective of the (poorly understood) conversion 
mechanism, the maximum energy is limited by the value 
7max = (e^/2meC^)iVGj,yiK^"^ (scc Lyubarsky & Kirk 
2001), which gives 7max = 0.9 x IQ^^k^^ and 1.6 x lOV;^^ 
for J1537A and J0737A, respectively. These limiting en- 



ergies are sufficiently large only at very small multiplici- 
ties (k ^ 1 for reasonable cr, n, and w), much lower that 
K ~ 100 expected from the pulsar models. Too low elec- 
tron energies and magnetic fields may be responsible for 
the fact that no TSs around supersonically moving MSPs 
have been convincingly detected in X-rays (e.g., Zavlin 
et al. 2002), in contrast to forward bow shocks seen in 

It is worth mentioning another problem with the wind- 
ISM interaction as the source of X-ray emission from 
close binary pulsars, such as J1537 and J0737. The con- 
ventional model for the shock (and the above consider- 
ation) assume that the pulsar's velocity relative to ISM 
is constant. However, a binary pulsar rotates around the 
center of mass, and the velocity of this motion may even 
exceed the binary systemic velocity. For instance, the 
A's orbital velocity varies between 150 and 263 km 
in J1537 and between 275 and 328 km s^^ for J0737. 
This means that both the absolute magnitude and the 
direction of the A's velocity relative to ISM vary with 
the orbital period. Moreover, using the pulsar's velocity 
relative to ISM in the above consideration would result 
in TS changing its position and orientation in the ISM 
(e.g., the TS head could appear behind the moving binary 
during some fraction of the binary period) . On the other 
hand, the characteristic time of propagation of pertur- 
bations in the ultrarelativistic wind, > R\i/c ^ 30-100 
hours, exceeds the binary period, so that we do not ex- 
pect the wind-ISM interaction to create a rotating (or 
wobbling) TS "bullet". Nevertheless, even if this inter- 
action results in a TS whose global properties, averaged 
over the orbital period, resemble those of a solitary pul- 
sar moving with a constant speed w, this structure should 
respond to the perturbations caused by the fast binary 
motion, and these perturbations should manifest them- 
selves in periodic variations in the radiation from the 
shocked plasma. However, no siich variations have been 
seen in the X-ray radiation from J0737. If future X-ray 
observations of J0737 put more stringent limit on the or- 
bital variation, it will be an additional argument against 
the interpretation in terms of the wind-ISM interaction. 

5. CONCLUSIONS 

We have detected the DNSB J1537 in X-rays, analyzed 
its X-ray spectrum and light curve, and compared the 
properties of its X-ray emission with those of the double 
pulsar J0737. Although the X-ray spectra and luminosi- 
ties of these systems are similar, J1537 shows a gap in 
its light curve during an interval of 0.35 of orbital period 
around apastron while J0737 does not show show a sig- 
nificant orbital phase dependence. There are two viable 
interpretations of the X-ray emission: it can be caused by 
interaction of the pulsar A's wind with pulsar B (§4.1), 
including the capture of the wind particles (§4.1.2), or it 
can be emission from the magnetosphere or, more likely, 
hot polar caps of pulsar A ( similar to solitary MSPs 
with E < lO"*^ ergs s"'^; §4.2). The former interpreta- 
tion can explain the putative orbital phase dependence 
in J1537 if the A's wind is mostly concentrated in the 
equatorial plane around the A's spin axis inclined to the 
orbital plane. The phase dependence can also be asso- 
ciated with the large eccentricity of J1537 (much larger 
than in J0737) because the X-ray luminosity depends on 
binary separation. On the contrary, no orbital phase de- 
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pendence is expected if the X-rays are produced by pulsar 
A. 

Unfortunately, the scarce count statistics does not al- 
low us to discriminate between the different interpreta- 
tions. Further observations of both J1537 and J0737, 
which have different orbital parameters and orientations 
with respect to the observer, should provide a clue to 
the nature of the X-ray radiation from DNSBs and the 
properties of pulsars and their winds. 

If future observations of J1537 confirm the deficit of 
X-ray flux during a substantial part of orbit, this would 
prove that the X-rays are produced via interaction be- 
tween the binary companions. If a similar deficit is ob- 
served at the opposite part of the orbit, it would strongly 
support the hypothesis that the A's wind is concentrated 
towards the equatorial plane, and it would allow one to 
measure the opening angle of the equatorial outflow. If, 
on the contrary, the X-ray flux shows a maximum at pe- 
riastron, it would mean that the binary separation is the 
main cause of the variable X-ray flux. For J0737, the 
dependence of the X-ray flux on orbital phase should be 
much weaker due to a smaller angle between the orbital 
angular momentum and the A's spin or/and the much 
smaller eccentricity. Confirming the interaction with the 
B's outer magnetosphere as the source of X-rays would 
provide strong evidence that the A's wind is accelerated 
up to energies > lO^meC^ at the interaction site, likely 
by magnetic reconnection. In such a scenario, one can 
expect X-ray pulsations with B's period, especially if the 
wind is captured by pulsar B. If the emission spectrum 
is found to be softer at the pulsations maxima, we would 
conclude that a fraction of the captured A's wind is chan- 
neled onto B's polar caps, heating them to X-ray tem- 
peratures. Detection of radiation caused by the capture 
of the A's wind would put a stringent lower limit on pair 
multiplicities in pulsar A. 

If future observations rule out any orbital phase de- 



pendence in the DNSBs (including the flux deficit near 
apastron in J1537), but pulsations at the A's period are 
detected, then we will have to conclude that the X-ray 
emission is produced by pulsar A itself, while the interac- 
tion of the A's wind with pulsar B is insignificant. This 
would mean that conversion of the intercepted wind's 
power into X-ray luminosity is inefficient, possibly be- 
cause the particles are not accelerated to high enough 
energies (which results in small values of ex and Crad in 
eq. [8]). In this case, studying the X-ray spectrum and 
pulsations would better characterize pulsar A, allowing 
one to compare its properties with those of "ordinary" 
MSPs, which have shorter periods and larger character- 
istic ages. 

Finally, it is quite plausible that the observed emission 
consists of two components: a soft thermal component 
from pulsar A and a hard nonthermal component from 
the A's wind interaction with pulsar B, as suggested by 
the X-ray spectra of J0737 and J1537. If this is the case, 
we expect to see periodic variations with both the or- 
bital and A's periods. Also, the relative contributions of 
these two components are expected to vary with orbital 
phase. The changes can be seen through phase-resolved 
spectroscopy because the wind emission is expected to 
be substantially harder than the MSP emission. 

To conclude, further deep observations of J1537 and 
J0737, and perhaps other DNSBs, would be of great in- 
terest as they will constrain both the properties of pulsar 
winds, including their interaction with the neutron star 
companions, and the properties of pulsars themselves. 
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TABLE 1 

Observed and derived parameters for J1537+1155 and J0737-3039 DNSBs 



Parameter 

Binary period, (days) 

Relative semimajor axis, a (10^^ cm) 

Eccentricity, e 

Distance, d (kpc) 

Orbital inclination, sini 

Transverse velocity, v±_ (km/s) 

Mass, M (Mq) 



J1537+1155 



J0737-3039 



0.421 

2.28 

0.274 

1.0 

0.975 

120 

2.68 



0.102 

0.878 

0.0877 

« 0.5 

1.000 

70 

2.59 



Note. — Based on the data from Stairs et al. (2002, 2004), Konacki et al. (2003), Lyne et al. (2004), Manchester et al. (2005), and Cole et al. 
(2005). The transverse velocoty in the plane of the sky, v±, is uncorrected for the solar system motion. 



TABLE 2 

Observed and derived parameters for J1537+1155A, J0737-3039A and J0737-3039B pulsars. 



Parameter 


J1537+1155A 


J0737-3039A 


J0737-3039B 


Spin period, P (ms) 


37.90 


22.70 


2773.46 


Period derivative, P (10"^*^) 


0.242 


0.174 


88 


Dispersion measure, DM (cm"* pc) ... 


11.6 


48.9 


48.7 


Angle S (deg) 


25 ± 4 (or 155 ± 4) 


0-60 (or 120 - 180) 




Angle a (deg) 


103 


0-60 




Mass (Mq) 


1.33 


1.34 


1.25 


Semimajor axis, a (10^^ cm) 


1.12 


0.424 


0.454 


Surface magnetic field, Bs (10^" G).... 


0.97 


0.64 


158 


Spin-down power, E (10** erg s~^).... 


1.8 


5.8 


0.0016 


Age, r - P/{2P), (Myr) 


248 


210 


50 


Light cylinder radius, Ri^. (10* cm) ... 


1.8 


1.1 


130 


Magnetic field at Ric, Bic (G) 


1660 


5050 


0.69 


Polar cap radius, Rpc (m) 


740 


960 


87 


Primary pair flux, Nqj (10**^ s~^) 


19 


34 


0.6 



Note. — Based on the data from Stairs et al. (2002, 2004), Lyne et al. (2004), and Manchester et al. (2005). Angle S is between the pulsar's 
spin and orbital rotation axis. Angle a is between the pulsar's spin and magnetic axis. The surface magnetic field, Bs, is assumed to be that of a 
dipole and given at the equator. 



TABLE 3 

Parameters of PL and BB models fitted to the Chandra and XMM-Newton spectra of J1537 and J0737 



Instrument Model Nh,2o'^ 


Norm.'' 


F or kT" 


C/dof 




Lx or Lbol(cos?)'' 


ACIS PL 3.6 
ACIS BB 3.6 


J1537+1155 

5.7li;^ X lO-'' 
4.8l3;o X lO'^ 


, -17+0.52 
•^•^'-0.52 


1.1/2 
5.3/2 


12% 
76% 


6.1+1° X 10^9 
e.Stil X 10^8 



ACIS 


PL 


5.0 


i.oii: 


ACIS 


BB 


5.0 


1.9lj 


MOSl+2 


PL 


5.0 


0.9811 


MOSl+2 


BB 


5.0 


3.95' 



0.16 
fO.99 



J0737-3039 

" ' X lO"-""' 
X 10* 
X 10-5 
X 10* 



0.08 
0.10 
+1.5 
1.0 



0.29 
0.31 
-17 
-15 
f 0.19 
-0.19 

16911^ 



2.93: 
192 
3.35' 



1.7/4 5% 



+0.42 
-0.39 

5.7/4 52% 2.67t.o lr ^ 
10.0/14 12% 3.03lo;49 X 10* 
16.9/14 62% 3.28lg j° x 10^ 



Note. — The errors are at the 68% confidence level for one interesting parameter. 
^Fixed values of hydrogen column density, in units of 10^*^ cm"^. 

^Normalization: Spectral flux in photons cm" ^ s" ^ keV" ^ at 1 keV or projected area of the emitting region in cm^, for the PL and BB models, 
respectively. 
^Photon index or BB temperature in eV. 

*^PercGntage of 10,000 Monte Carlo simulations, drawn from the best-fit model, whicli give a C-statistic value lower than the best-fit value. Although 
the BB fits arc formally worse than the PL fits, they are statistically acceptable. 
^ Lx is the unabsorbed PL luminosity in the 0.2—10 keV band, -Lboi the bolometric luminosity for the BB fits. 



